Quantum Zeno Effect in the Quantum Non-Demolition Detection of Itinerant Photons 



00 
O 

o 

(N 



O 



a 

•4-* 

s 

C 

o 
o 



(N 
> 

oo 
o 

OS 

o 



X 



Ferdinand Helmer, 1 Matteo Mariantoni, 2 Enrique Solano, 1 and Florian Marquardt 1 

Department of Physics, Center for NanoScience, 
and Arnold- Sommerf eld Center for Theoretical Physics, 
Ludwig-Maxirnilians-Universitdt, Theresienstrasse 37, D-80333 Munich, German^ 
2 Walther-Meissner-Institut, Bayerische Akaderaie der Wissenschaften, 
Walther-Meissner-Strasse 8, D-85748 Garching, Germany 

We analyze the detection of itinerant photons using a quantum non-demolition (QND) measure- 
ment. We show that the backaction due to the continuous measurement imposes a limit on the 
detector efficiency in such a scheme. We illustrate this using a setup where signal photons have to 
enter a cavity in order to be detected dispersively. In this approach, the measurement signal is the 
phase shift imparted to an intense beam passing through a second cavity mode. The restrictions on 
the fidelity are a consequence of the Quantum Zeno effect, and we discuss both analytical results 
and quantum trajectory simulations of the measurement process. 



Quantum non-demolition (QND) measurements are 
ideal projective measurements that reproduce their out- 
come when repeated 0,0. Using them, it is possible to 
measure the state of a system with the minimal disrup- 
tion required by quantum mechanics. Recent successful 
experimental demonstrations of QND detection for su- 
perconducting qubits and microwave photons [1, 0> H @| 
are both of fundamental interest and crucial for the de- 
velopment of quantum communication and information 
processing. When QND detection is applied continuously 
to a system that would otherwise undergo some intrin- 
sic dynamics, quantum jumps are observed, tracing the 
quantum evolution in real-time 0, S, HI EH EH ■ As a con- 
sequence, the dynamics tends to be frozen, a result now 
lHj I known as the Quantum Zeno effect [H, El El EH 0, El | • 



In the present paper, we show that the interplay of 
these phenomena may put interesting constraints on the 
detection of itinerant quanta. The specific minimal ex- 
ample we will discuss concerns the continuous dispersive 
QND detection of single photons passing through a cav- 
ity. The crucial distinction to be recognized is the follow- 
ing: For localized quanta (e.g. a photon already created 
inside a cavity [H EH), the Quantum Zeno effect could 
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Figure 1: (Color online) Schematic sketch of the model. Two 
cavity modes are coupled anharmonically. The detector mode 
is irradiated with a strong coherent field that suffers a phase 
shift whenever a photon is present in the signal mode. 



presumably only enhance the detection by suppressing 
the decay. However, this no longer holds for the detec- 
tion of itinerant quanta, if we require that our detector 
is always working and can detect the quantum without 
knowing the arrival time in advance. As we will show, in 
this case the unavoidable back-action of the measurement 
device produces a Quantum Zeno effect, suppressing the 
fidelity of measurements. 

We investigate a QND scheme utilizing the non- linear, 
2(il . 21 1 of two discrete localized 
The presence of a quantum 



Kerr-type coupling 1 



modes of a bosonic field, 
inside the signal mode gives rise to a frequency shift of 
the detection mode, which can be observed dispersively 
via the phase shift of a beam transmitted through that 
mode (see Fig. Q]). In turn, the signal mode frequency 
fluctuates due to the detection beam's shot noise. As a 
consequence, the incoming signal photon will be reflected 
with a probability that rises with coupling strength and 
detection beam intensity. 

This incarnation of the Quantum Zeno effect generates 
a trade-off that yields the highest detection efficiency at 
intermediate coupling strengths. In that way, such dis- 
persive schemes for itinerant quanta turn out to be simi- 
lar to weak measurements using general linear detectors 
and amplifiers Q. 

We proceed as follows: (i) We numerically evaluate 
quantum jump trajectories for the phase shift signal in 
a minimal model of a QND photon detector and analyze 
the fraction of detected photons, observing the trade-off 
described above, (ii) We interpret these findings using an 
analytical approximation, (iii) Finally, we briefly com- 
ment on possible experimental realizations. 

Model. - We consider a system of two cavity modes 
with a Kerr-type coupling of strength g: 

H = huj 
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kg h TT-dct + -ffdrivc+dccay ■ (1) 

These modes might represent two different electromag- 
netic field modes inside an optical or microwave cavity, 



2 




10 20 30 

time K,t 



40 50 



Figure 2: (Color online) Quantum jump trajectories illustrat- 
ing dispersive photon detection: The observable homodyne 
signal X(t) (red lines) and the corresponding signal mode oc- 
cupation (n) (t) (blue lines), for two different values of the 
measurement rate T/n at a fixed input rate Ni n - Photon de- 
tection events are indicated as filled circles. The relative noise 
strength oc (rr aV g)~ 1//2 is suppressed with increasing T/k, but 
the number of photons actually detected also decreases, due 
to the Quantum Zeno effect (see main text). The size of the 
noise floor, the detector threshold X t h r (green dashed line), 
and the dark time are indicated. Here and in the following 

plots KT avg = 2. 



the modes of two adjacent cavities, or even two anhar- 
monically coupled modes of a nanomechanical resonator. 
Photons in the signal mode (frequency u), number op- 
erator n) and the detector mode (wdct, ^det) decay by 
leaking out of the cavity. The anharmonic Kerr-type cou- 
pling arises generically when introducing any nonlinear 
medium, such as an atom, a qubit or a quantum dot, into 
a cavity and has been studied for the purpose of QND 
measurements in quantum optics H, 2(J, HH . It induces 
a phase shift in the strong detection beam ({hdet) ^ 1) 
upon presence of a signal photon. 

We are interested in analyzing individual realizations 
of the phase shift signal as a function of time. The phase 
shift can be observed by continuously measuring an ap- 
propriate field quadrature of the detection beam (e.g. in 
a homodyne setup) . As the beam passes through the cav- 
ity, the beam becomes entangled weakly with the cavity's 
state. Thus, the stochastic measurement outcomes reveal 
information about that state, feeding back into the time- 
evolution of the cavity's density matrix. This physics 



is described by a stochastic Lindblad master equation 
@,[H, H, 23|, 24 1 for the density matrix p conditioned on 
the output signal (see fl8j|): 
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-2r [n, [h, p]] - Vlf (hp + ph-2p (h) (t)) f (t).(2) 

We analyze a situation with a continuous weak coherent 
beam of photons entering at a rate N- m into the signal 
mode, whose decay-rate is k (first line of Eq. ((2])). We 
have chosen to work in the limit of a large detector mode 
decay rate, Kdot 3> K, which is favorable for the detection 
process and makes it possible to adiabatically eliminate 



that mode [la ], keeping only the signal mode h = tva and 



drastically reducing the numerical effort. After adiabatic 
elimination, the coupling strength g and the detection 
beam intensity are combined into the measurement rate 
18] r = g 2 (ndct) /(4«dct), where 1/r is the time-scale 
needed to resolve different photon numbers. The last, 
stochastic term in Eq. describes the measurement 
back action. The (suitably normalized) phase shift signal 
reads: 



(3) 



It contains a systematic term depending on the average 
number of signal photons, as well as a stochastic term 
representing the unavoidable vacuum noise, where (£) = 
and (£(*)£(*')) = 8{t - t'). In deriving Eq. ©, we 
have assumed that the transmitted and reflected signals 
are superimposed symmetrically to extract the maximum 
information. As in any measurement of field quadratures, 
temporal filtering is required to suppress the noise. We 
average over a timespan T avg , which should be as large as 
possible while still remaining smaller than the expected 
temporal extent of the phase shift signal due to a single 
photon, i.e. r avg <C k^ 1 . We denote the averaged signal 
as X(t). 

Numerical Results. - We numerically solve the master 
equation, using it to compute the signal X(t) and the oc- 
cupation of the signal mode (h) (t) as a function of time. 
We then implement the minimal model of a threshold 
detector: Time-points when the quantum jump trajec- 
tory X(t) first exceeds the threshold A t h r are counted as 
detection events, and the detector is then set insensitive 
for a dark-time Tdarkj suitably chosen to avoid multiple 
detection, i.e. N^ 1 » ra a rk ^ Our discussion will 

focus on small values of N[ n , making the results inde- 
pendent of Tdark, while we will analyze the dependence 
on A t hr in some detail. In figure [2] we show two exam- 
ple trajectories. Whereas the expected number of signal 
photons is the same for both cases, the increase in the 
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measurement rate Y/k decreases the number of photons 
actually detected, while at the same time enhancing the 
signal to noise ratio in the trajectory. This is an indi- 
cation of the Quantum Zeno effect in the detection of 
itinerant quanta, which we now want to study in more 
quantitative detail. 

We plot the rate of photon detection events Ndet ver- 
sus the rate of incoming photons N- m [Fig. OH a)]. The 
detection efficiency rj is naturally defined as the ratio of 
detected vs. incoming photons, obtained at small input 
rates iV; n : 



dK 



det 



(4) 
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Fig-0H c ) displays the efficiency 77 as a function of Y/k and 
X-thi- The statistics for this figure were obtained from ex- 
tensive numerical simulations by generating O(10 4 ) tra- 
jectories of length 10 2 /k for seven different rates N- m at 
each value of Y/k. Apparently, the detector efficiency 77 
is strongly suppressed both for Y/k <C 1 (low signal to 
noise ratio) and Y/k ^> f . 

Analytical Results. - To interpret these results, we now 
calculate the total transmission probability through the 
signal mode, whose frequency fluctuates due to the shot 
noise in the detection mode, which is treated as classical 
noise. We start from the semiclassical equation of motion 
for the complex field amplitude a(t) in the signal mode, 



d(t) 



-{ du)(t) 



a(t) 



(5) 



Here a™ is the amplitude of the signal photon field en- 
tering the cavity from the left side, and 5oj(t) = gridet(t) 
is the fluctuating frequency shift (ridct S> 1). The corre- 
lator of the noise is given by 



(6) 



To obtain an expression for the transmission probability, 
we write down the formal solution for a(t), 



fKLCt]' 



dt! exp 



Su{t")dt" -~{t 







(7) 

Note that the fluctuations Su>(t) themselves are non- 
Gaussian. Still, the integral in the exponent is approx- 
imately Gaussian for time- intervals that fulfill Kdct|£ — 
t'\ ^> f, due to the central limit theorem. These times 
yield the main contribution under our assumption of a 
"fast detector", K^ct ^> K. Thus, we can evaluate (|ev| 2 ) 
using the formula (exp[— iY]) = exp[—i(Y) — ^VarF] for 
a Gaussian random variable Y, and inserting Eq. ([6]). 
From this, we obtain the average transmitted intensity 
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Figure 3: (Color online) (a) Detector profile: Rate of detected 
vs. incoming photons, at T/k = 0.6, for two different thresh- 
olds Xthi- Observe the dark count rate (offset at jVj n = 0), the 
detector efficiency n defined from the slope at N- m = 0, and 
the saturation for large N- ln ~ ^dark- ( D ) Suppression of the 
signal photon number inside the cavity as a function of mea- 
surement rate Y/k. (c) Detector efficiency n, obtained from 
quantum trajectory simulations, as a function of Y and X t hr 
(3D-inset), and comparison with the analytical results (main 
plot). The blue data points display r\ for fixed — 0.5 

(blue cut in 3D-inset). When maximizing the efficiency over 
X t hr for any given Y/k, the red data points are obtained (la- 
beled "max(?7)"), coinciding with the analytical asymptote 
(green thin line) at higher values of Y/k. Small inset: The 
dependence of the efficiency n on the averaging time r avg is 
shown for Y/k = 4. 



and the average transmission probability 
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Before we can correlate the suppression of the transmis- 
sion with the reduction of the detector efficiency rj in the 
limit of T/k > f , one more consideration is necessary. 



4 



In this limit, any photon that has entered the cavity will 
almost certainly be detected. Once detected, the pho- 
ton loses the coherence with the incoming beam, which 
is needed for perfect transmission on resonance in the 
ideal, coherent case. As a consequence, it acquires an 
equal probability to leave the cavity through the left or 
the right port. This means that, on average, the number 
of detected photons is twice the number of transmitted 
photons. The expected relation is thus n = 2(T), which 
is indeed observed nicely when comparing with the nu- 
merical data (Fig. EUc)). The reduction of detector ef- 
ficiency at T/k > 1 thus has found its explanation in 
the Quantum Zeno effect: many photons remain unde- 
tected, because they are reflected due to detector back- 
action. As low values of T/k are also unfavorable, due 
to a bad signal-to-noise ratio, the maximum efficiency is 
found at an intermediate value, namely near T/k = 1/2 
(see Fig. |3J}. 

Possible realizations. - Cavity QED setups in super- 
conducting circuits [12, 25, [3, [53] have been used to 
implement ideas of quantum optics on the chip, and are 
considered a promising candidate for scalable, fault tol- 
erant quantum computing While proposals for gen- 
erating nonclassical photon states exist or have been im- 
plemented 0, 0, 23, 3(| , the on-chip single-shot detec- 
tion of itinerant photons is still missing. Building on 
recent experiments that demonstrated dispersive qubit 
detection Q and measurements of photon statistics 
one could employ the superconducting qubit to induce a 
nonlinear coupling between two modes of the microwave 
transmission line resonator (or coupling two cavities [33j]). 
thus creating a dispersive photon detector of the type 
discussed here. These experiments realize a Jaynes- 
Cummings coupling between qubit and resonator of up 
to 2ir ■ 100MHz, resonators with frequencies of about 
2ir ■ 5GHz, and a large spread of resonator decay rates 
k between 1MHz and 100MHz. Given this wide param- 
eter range, it is possible to cover the full range of T/k 
explored here, assuming around 10 to 100 photons in the 
detection mode. The detector efficiency, although limited 
by the Quantum Zeno effect as shown before, could then 
reach values of about 40% even without considering more 
elaborate detector and signal analysis schemes. Another 
experiment in which essentially the same physics could 
be observed is the detection of single photons in a mi- 
crowave cavity by employing the dispersive interaction 
with a stream of Rydberg atoms Q . 

Conclusions. - In this paper we have analyzed a rather 
generic scheme for the detection of itinerant photons in 
a QND measurement process, employing quantum tra- 
jectory simulations. We have shown how the Quantum 
Zeno effect enters the detection efficiency, a result that 
will be relevant to many other situations, such as the 
detection of electrons tunneling through a quantum dot 
by current passing through a nearby quantum point con- 
tact 3l|, the detection of itinerant phonons entering a 



micromechanical cantilever or membrane (e.g. in an op- 
tomechanical setup [lH), and other similar settings in 
mesoscopic physics, quantum optics, and atomic physics. 
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